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Resumen 


Larios-Tlali, H., Torres-Benites, E., Quevedo-Nolascol, A., 
Martínez-Menes, M. R., € Salgado-Tránsito, J. A. (mayo- 
junio, 2015). Riesgo de inundación en la subcuenca del río 
La Antigua, Veracruz, México. Tecnología y Ciencias del Agua, 
6(3), 39-56. 


El objetivo del presente trabajo fue aplicar el modelo de 
simulación hidráulica HEC-RAS para delimitar las zonas 
con riesgo a inundación en la subcuenca del río La Antigua, 
Veracruz, México, utilizando la metodología del Centro 
Nacional de Prevención de Desastres (Cenapred), que 
considera la vulnerabilidad y el peligro de un sistema ante la 
ocurrencia de un fenómeno extremo perturbador. El modelo 
HEC-RAS usa ecuaciones diferenciales determinísticas que 
simulan el comportamiento y la dinámica de los niveles 
de agua en secciones transversales del cauce; su aplicación 
bajo condiciones de flujo no permanente y régimen mixto 
permitió definir las áreas de inundación alcanzadas 
durante una avenida extraordinaria. Para la simulación 
hidráulica se construyeron hidrogramas unitarios sintéticos 
con datos de precipitación registrados durante un evento 
hidrometeorológico extremo; la geometría del cauce se 
obtuvo con la extensión HEC-GeoRAS, herramienta que 
extrajo y procesó información georreferenciada del Modelo 
Digital de Elevación (MDE) LiDAR tipo terreno y los 
coeficientes de rugosidad (n de Manning) se estimaron de 
acuerdo con las características físicas del canal. Para calibrar 
los datos estimados por el modelo se utilizó información 
hidrométrica (niveles registrados) de la estación 28003 
Cardel. La aplicación de tres estadísticas de prueba: R, RMSE 
y DMA demostraron un alto ajuste entre datos observados 
versus simulados, respaldando la aplicabilidad del modelo 
HEC-RAS para simular, predecir y delimitar zonas de riesgo 
a inundación. 

Palabras clave: simulación hidráulica, inundación, 
hidrograma, fenómenos hidrometeorológicos extremos, 
calibración, riesgo a inundación. 


Abstract 


Larios-Tlali, H., Torres-Benites, E., Quevedo-Nolascol, A., 
Martínez-Menes, M. R., €' Salgado-Tránsito, J. A. (May-June, 
2015). Flood Risk in a Sub-Basin of the La Antigua River, Veracruz. 
Mexico. Water Technology and Sciences (in Spanish), 6(3), 39- 
56. 


The aim of this study was to apply the HEC-RAS hydraulic 
simulation model to delineate flood hazard areas in a sub-basin of the 
La Antigua River, Veracruz, using the methodology by the National 
Center for Disaster Prevention (Cenapred). This methodology takes 
into account the vulnerability and risk of a system to an extreme 
disturbance. The HEC-RAS model uses deterministic differential 
equations to simulate the behavior and dynamics of water levels 
in cross-sections of channels. By applying it under unsteady and 
mixed flow conditions, areas of flooding resulting from an unusually 
large flood can be defined. For hydraulic simulation, synthetic 
unit hydrographs were constructed using data recorded during 
an extreme hydrometeorological precipitation event. Channel 
geometry was obtained with HEC-GeoRAS, a tool which extracts 
and processes georeferenced information from a LiDAR land 
surface Digital Elevation Model (DEM). Roughness coefficients (n 
Manning) were estimated according to the physical characteristics 
of the channel. To calibrate the estimates calculated by the model, 
hydrometric information (recorded levels) was used from the Cardel 
28003 station. The application of three test statistics— R, RMSE 
and DMA— showed a good fit between simulated and observed 
data, indicating the applicability of the HEC-RAS model to simulate, 
predict and delineate flood hazard zones. 


Keywords: Hydraulic simulation, flood hydrograph, extreme 
hydrometeorological events, calibration, flood. 
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Introducción 


México, por su ubicación geográfica, es un país 
vulnerable a la ocurrencia de diversos fenóme- 
nos hidrometeorológicos, principalmente a la 
formación y desarrollo de ciclones tropicales 
(Jáuregui, 2003) en la costa del Pacífico y del 
Golfo de México. Un efecto de estos fenómenos 
—en cualquiera de sus etapas— son las inunda- 
ciones originadas por lluvias torrenciales; sin 
embargo, no es el único factor que determina 
la magnitud de estas catástrofes, también 
incluye la distribución espacial de la lluvia, la 
morfometría de los cauces, la infraestructura, la 
topografía, el relieve, la cobertura vegetal, el uso 
de suelo y las unidades edafológicas (Cenapred, 
2001). 

La zona costera del estado de Veracruz, 
por su ubicación geográfica, orografía, clima 
y asentamientos humanos en zonas de alto 
riesgo, también es un área vulnerable a estos 
fenómenos. En las dos últimas décadas (de 1990 
a 2010), los municipios de Tuxpan, Tamiahua, 
Tecolutla, Nautla, San Andrés Tuxtla y La An- 
tigua sufrieron los efectos de cinco tormentas 
tropicales y diez huracanes (GEV, 2011). 

El paso cercano de tormentas o ciclones 
tropicales pueden ser causas de inundaciones, 
pero también lo pueden ser las depresiones 
tropicales menores e incluso nortes (Acevedo éz 
Luna-Díaz, 2006). Entre 1950 y 1994 se registra- 
ron 475 inundaciones en la región Golfo Centro, 
de las cuales 255 ocurrieron en la región Pánuco 
(CSVA, 2007), la mayoría de estas catástrofes 
se relacionan con la ocurrencia de fenómenos 
hidrometeorológicos extremos. 

El potencial hídrico superficial del estado 
de Veracruz es de los más altos de México, su 
escurrimiento superficial anual medio se esti- 
ma en 121 millones de metros cúbicos, lo que 
representa 33% del escurrimiento superficial de 
todo el país (CSVA, 2005). Veracruz tiene una su- 
perficie de 72 815 km?, 73% de su territorio está 
por debajo de los 200 msnm (Tejeda, 2006). El 
gobierno del estado de Veracruz (GEV, 2011) es- 
tima una superficie inundable de 6 275 km? (8% 
del territorio estatal), donde habitan 1 252 251 


personas en 2 189 localidades de 118 municipios, 
lo que indica que uno de cada seis veracruzanos 
se ve afectado por las inundaciones. 

Por lo anterior, es necesario e imprescindible 
realizar un estudio integral para estas zonas con 
alto riesgo a inundación, el cual debe proponer 
acciones, Obras y programas que prevengan, 
reduzcan, atiendan y mitiguen los daños ocasio- 
nados por tales catástrofes. Una de las medidas 
no estructurales para disminuir el riesgo por 
inundación es la delimitación de dichas zonas. 

El presente trabajo ha combinado el uso 
de sistemas de información geográfica (SIG), 
técnicas de percepción remota, trabajos de 
prospección en campo, análisis de información 
hidrometeorológica y aplicación de un modelo 
matemático de simulación hidráulica (HEC- 
RAS, versión 4.1.0) para delimitar las zonas 
con riesgo a inundación en la subcuenca del río 
La Antigua, Veracruz, una de las regiones del 
estado más dañadas por inundaciones en las 
últimos años. 

Los resultados obtenidos en la simulación 
hidráulica fueron comparados con los niveles 
registrados en la estación hidrométrica 28003 
Cardel (durante la ocurrencia de un fenómeno 
hidrometeorológico extremo), para ello se em- 
plearon tres estadísticas de prueba (R, RMSE y 
DMA), cuyos valores obtenidos demostraron 
un alto ajuste entre datos observados versus 
simulados. 


Área de estudio 


La cuenca del río La Antigua se ubica entre los 
paralelos 19? 05' y 19? 34” latitud norte y los 
meridianos -96* 06' y -97* 16' longitud oeste. 
Tiene un área aproximada de 2 827 km?, que se 
distribuye en los estados de Puebla y Veracruz. 
El río La Antigua nace en la Sierra Madre Orien- 
tal, con el nombre de río Sumidero, al oriente 
de la localidad González Ortega, en el estado 
de Puebla, a una altitud de 3 345 msnm (Espi- 
noza, 2009). La cuenca pertenece a la Región 
Hidrológica 28 Papaloapan (RH28) y a la Su- 
bregión Jamapa y Otros, y está conformada por 
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cuatro subcuencas: RH28Bc, RH28Bd, RH28Be 
y RH28Bf (figura 1). 

El área de estudio incluye sólo la subcuenca 
RH28Be (subcuenca del río La Antigua), ya que 
en esta zona se registran los mayores daños 
provocados por el desbordamiento del río La 
Antigua. 

La subcuenca RH28Be tiene una superficie 
de 508.3 km?, se enmarca dentro de los paralelos 
19* 12 y 19* 23' y los meridianos -96* 21” y -97* 
00”. Se distribuye horizontalmente con dirección 
oeste-este en el estado de Veracruz, sobre las 
subprovincias fisiográficas denominadas “Chi- 
conquico” y “Llanura Costera Veracruzana”, las 
cuales pertenecen, respectivamente, a la Provin- 
cia Fisiográfica X Eje Neovolcánico Transversal 
y a la Provincia Fisiográfica XIII Llanura Costera 
del Golfo Sur. 

El clima predominante, según la clasificación 
de Kóppen, modificada por Enriqueta García, 


es el grupo A (cálidos húmedos tropicales), 
donde domina la unidad climática “Awo”, con 
una temperatura que oscila entre 18.7” y 25.2 *C; 
en los meses de abril a mayo se presentan los 
valores más altos y de noviembre a febrero las 
temperaturas más frías; la temperatura media 
anual es de 22.3 *C. Por otro lado, la precipita- 
ción media anual es de 1 393 mm, en la parte 
alta de la subcuenca (zona montañosa) se han 
registrado valores que superan los 2 000 mm, y 
en la parte central se registran valores entre 900 
y 1000 mm. La evaporación media anual es de 
1 340 mm, en los meses de marzo a mayo es más 
intensa y en el mes de diciembre se presentan 
los valores mínimos (SMN, 2008). De acuerdo 
con el conjunto de datos vectoriales de uso de 
suelo y vegetación serie IV, escala 1:250 000 del 
INEGI (2012), en la subcuenca RH28Be existen 
seis usos de suelo y vegetación: agricultura, 
pastizal cultivado, bosque de encino, bosque de 
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Figura 1. Ubicación geográfica de la cuenca del río La Antigua, con base en la Red Hidrográfica escala 1:50 000 del Instituto 
Nacional de Estadística y Geografía (INEGD), 2010. 
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coníferas, selva baja caducifolia y selva mediana 
subcaducifolia; el uso de suelo predominante 
es la agricultura, ocupa una superficie de 322.3 
km, lo que representa 63.4% del área total. Por 
otro lado, con base en el Sistema de Clasificación 
de Suelos WRB-2000, la unidad edafológica 
de mayor distribución es el vertisol pélico, al 
ocupar una superficie de 271.7 km”, lo que re- 
presenta 53%, seguido del acrisol órtico y litosol, 
con 110.7 y 98.3 km”, respectivamente. 


Materiales y métodos 
Descripción del modelo 


El HEC-RAS (Hydrologic Engineering Center- 
River Analisis System) es un modelo numé- 
rico creado por el Cuerpo de Ingenieros de la 
Armada de los Estados Unidos, de uso libre 
y en continuo desarrollo (Turégano, 2012). El 
procedimiento básico de cálculo de este software 
se fundamenta en la solución de la ecuación de 
energía, donde las pérdidas de energía por fric- 
ción se estiman usando la ecuación de Manning; 
las pérdidas locales por contracción y expansión 
del flujo se calculan como una fracción del cam- 
bio en la velocidad entre dos secciones (USACE, 
2002). Este programa utiliza el método de paso 
estándar (Standard Step Method), el cual consis- 
te en determinar el tirante en una sección a partir 
del tirante asociado con el escurrimiento de una 
sección adyacente (aguas arriba o aguas abajo, 
dependiendo del régimen de escurrimiento). El 
balance de energía entre estas dos secciones se 
hace resolviendo la siguiente ecuación: 


donde Y,, Y, = tirante en la sección trasversal; 
Z,, Z, = elevaciones del fondo del canal prin- 
cipal; V,, V, = velocidades promedio; a, 4, = 
coeficientes de velocidad; g = aceleración de la 
gravedad, y h, = pérdida de energía (Ceballos, 
2011). 

Según Timbadiya, Patel y Porey (2011), este 
razonamiento simplifica el flujo de agua como 


un flujo unidimensional, asume que el tirante y 
velocidad solo varían en dirección longitudinal 
del canal, cuyo eje se supone aproximadamente 
a una línea recta con velocidad constante en 
cualquier punto de la sección trasversal. 


Delimitación de las unidades de análisis 


La delimitación de la unidad de análisis esta- 
blece las interacciones entre los componentes 
del sistema y los procesos que en él ocurren 
durante la simulación hidráulica. El mapa base 
se elaboró utilizando los datos vectoriales de la 
RH28B, extraída de la Red Hidrológica escala 
1:50 000 reportada por INEGI en su página web 
(consultado el 1? de julio de 2013, disponible en 
http: / /www.inegi.org.mx). 

El análisis, como se ha señalado, incluyó sólo 
la subcuenca RH28Be, debido a que en esta zona 
se han registrado los mayores daños provocados 
por el desbordamiento del río La Antigua; sin 
embargo, considerando que el escurrimiento 
registrado por la estación hidrométrica 28003 es 
el acumulado de RH28Be y RH28Bf (figura 2), 
cuyos datos fueron utilizados en la calibración, 
se modeló el comportamiento y la dinámica de 
la red hidrológica RH28Be, tomando en cuenta 
el escurrimiento generado por RH28Bf. 

En la figura 3 se observa la suma del escu- 
rrimiento de RH28Bf al tramo “Sta. María”, a 
través del afluente conocido como “Pescados”. 
De esta forma se modelaron ocho cauces de 
RHA28Be y uno de RH28Bf (ver figura 2). 

Para cada afluente se requirió de informa- 
ción específica de su zona de dominio, por lo 
que se dividió al área de estudio en nueve uni- 
dades de respuesta hidrológica (URH); cada una 
corresponde a un afluente o tramo considerado 
en la modelación. Los parámetros de las URH 
que se utilizaron para estimar los caudales se 
muestran en el cuadro 1. 


Estimación de los hidrogramas de entrada 
En la subcuenca RH28Be, los registros de pre- 


cipitación son más abundantes que los de escu- 
rrimiento, por ello se eligió un método indirecto 
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Figura 2. Distribución de las subcuencas y unidades de respuesta hidrológica (URH). Para delimitar las URH se empleó WMS, 
ver. 8.3 y el MDE LiDAR tipo terreno del INEGI, 2013. 
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Figura 3. Confluencia del río Pescados (subcuenca RH28Bf) y el río La Antigua (subcuenca RH28Be). 


para construir un hidrograma de escorrentía fue el hidrograma adimensional curvilíneo del 
superficial, asociado con una lluvia de duración SCS (Soil Conservation Service), que puede ser 
y profundidad unitaria. El método utilizado representado por un hidrograma unitario trian- 
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Cuadro 1. Parámetros morfométricos de las URH. 


URH Parámetro 

A A car rar el al 
Xochilapa I 38.8 1556 260 193.96 0.032 
Actipan Tr 16.3 984 260 44.91 0.044 
La Bomba HI 35.8 900 182 53.96 0.020 
Huilapa IV 28.3 22 182 40.57 0.019 
Pescados v 113.9 3721 98 1 631.67 0.031 
Sta. María vI 63.0 1556 98 314.02 0.023 
La Bomba 2 vI 46.4 900 73 110.06 0.018 
Sta. María 2 VIH 120.0 3721 73) 1 950.06 0.030 
Parte Baja IX 142.2 3721 5 2 145.15 0.026 


gular (HUT). El método propuesto por Mockus 
(1957) genera un gasto pico en función del área 
de la cuenca y el tiempo pico según la siguiente 
ecuación: 


0.555 hp, A 
EN 


Pp 


2) 


Del análisis de varios hidrogramas, Mockus 
concluye que el tiempo base y el tiempo pico se 
relacionan mediante la expresión: 


t, =2.67 t, (3) 
el tiempo pico se expresa: 
d 
t, = a +1, (4) 
t,=0.6t (5) 


El tiempo de concentración (t) se puede 
estimar con la ecuación de Kirpich: 


0.77 
E =0.000325%— 


90.385 


(6) 


La duración del exceso de lluvia (d,) se 
calcula para cuencas grandes con la siguiente 
ecuación: 


(7) 


donde q, = gasto pico unitario (m*/s/mm); 


hp, = altura de precipitación efectiva (mm); 
A = área de la cuenca (km”); t, = tiempo pico 
(h); t, = tiempo base (h); t_ = tiempo retraso (h); 
L = longitud del cauce principal (m), y S = pen- 
diente del cauce (%.). 

Una vez generado el HUT con la metodolo- 
gía de Mockus, se estimó el hidrograma adimen- 
sional curvilíneo del SCS, que suaviza la forma 
del HUT: el hidrograma del SCS conserva las 
mismas unidades de tiempo y volumen que el 
HUT. Para obtener el hidrograma unitario del 
SCS se multiplican las ordenadas por el q, y las 
abscisas por el £, (Barrios $: Olaya, 2007). 

El hidrograma unitario se construyó con los 
datos de precipitación registrados en 25 esta- 
ciones meteorológicas (ubicadas dentro y en la 
periferia de la subcuenca RH28Be) durante la 
ocurrencia de un fenómeno hidrometeorológico 
extremo que impactó las costas de Veracruz el 17 
de septiembre 2010: el huracán Karl. Los datos 
de precipitación fueron proporcionados por la 
Dirección Técnica del Organismo de Cuenca 
Golfo Centro (OCGC) de la Comisión Nacional 
del Agua (Conagua). Con la información de las 
estaciones se calculó la precipitación media pon- 
derada de cada URH utilizando el método de las 
isoyetas, el cual establece un área de influencia 
delimitada por líneas que unen los puntos de 
igual valor de precipitación de acuerdo con la 
siguiente expresión: 


N 
D= gt (8) 
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donde D = altura de precipitación media (mm); 
a,= área entre cada dos isoyetas (km?); D,= pro- 
medio de precipitación entre dos isoyetas (mm), 
y A = área total de la cuenca (km). 

En la figura 4 se muestran los hidrogramas 
adimensionales curvilíneos del SCS estimados 
para las nueve URH, y en el cuadro 2 se tienen 
los valores de t, y Q, valores que construyen 
dichos hidrogramas. 


Geometría de cauces 


La geometría de los cauces se obtuvo modelan- 
do con HEC-GeoRAS, una extensión del HEC- 
RAS que extrae la información georreferenciada 
de un MDE. En este trabajo se utilizó el MDE 
LiDAR tipo terreno, con una resolución de 5 m 
en la horizontal y 0.25 m en la vertical. 

Se inició digitalizando el centro del cauce 
principal de la red hidrológica RHA28Be, comen- 
zando de aguas arriba hacia aguas abajo; los 
centros de cada afluente tributario se unieron 


conforme avanzó la digitalización. Posterior- 
mente se dibujaron las márgenes, primero 
la margen izquierda (del cauce principal) y 
después su margen derecha; de igual forma se 
unieron las márgenes de los afluentes tributarios 
de aguas arriba hacia agua abajo. La longitud de 
los centros de cauce digitalizados suman 186.5 
km y las márgenes suman 374.4 km (figura 5). 
Para generar las secciones transversales se 
tomaron en cuenta las siguientes condiciones: a) 
el trazo de la sección parte de la margen izquier- 
da hacia la margen derecha; b) las secciones se 
trazan de aguas arriba hacia aguas abajo y son 
perpendiculares al centro del cauce; c) las sec- 
ciones no se pueden interceptar entre ellas, y d) 
las secciones transversales sólo pueden cruzar 
una sola vez al centro y márgenes del cauce. En 
total se crearon 2 435 secciones transversales a lo 
largo de los nueve cauces o tramos que integran 
la red hidrológica digitalizada, el intervalo entre 
cada sección fue de 20 m y la longitud de trazo 
(o ancho de sección) osciló entre 100 y 1 200 m. 
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Figura 4. Hidrogramas adimensionales curvilíneos del SCS considerados en la simulación hidráulica. 
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Figura 5. Secciones transversales de los ríos La Bomba, Huilapa y La Bomba 2, obtenidas con la extensión HEC-GeoRAS. 


Determinación del coeficiente de Manning 


El coeficiente de rugosidad de Manning repre- 
senta la resistencia al flujo de agua en cauces 
sobre las llanuras de inundación; el valor es 
más alto cuando existe mayor rugosidad en la 
superficie de contacto con el flujo (Kumar, 2013). 
Este coeficiente debería ser calibrado a partir 
de información obtenida en campo (perfiles 
y secciones de ríos); sin embargo, cuando no 
se dispone de datos de aforo, se recomienda 
utilizar guías —para la selección del coeficien- 
te— que propongan valores para cauces de 
características similares o valores obtenidos de 
datos experimentales (Turégano, 2012). 

El modelo requiere la asignación de un coefi- 
ciente de rugosidad para la margen izquierda, 
derecha y centro de cada sección transversal. 
El valor se estableció utilizando la “Tabla del 
coeficiente de rugosidad de Mamning”, pro- 
puesta por Chow, Maidment y Mays (1994), la 
cual establece el coeficiente en relación con las 


características y condiciones más sobresalientes 
del canal, como rugosidad del fondo del cauce, 
cobertura y tipo de vegetación, y sinuosidad 
del canal, entre otras. La asignación de este 
coeficiente se apoyó en puntos de verificación 
en campo efectuados sobre el río La Antigua y 
sus afluentes. 

En la figura 6a se observa que el fondo del 
cauce posee distintas características de rugo- 
sidad; por lo tanto, el fondo de cada sección 
transversal posee un coeficiente de Manning 
acorde con las condiciones del sitio. También se 
observa que ambas márgenes pueden ser muy 
diferentes en una misma sección (figura 6b), por 
ello cada margen posee un coeficiente propio. 
En total se obtuvieron 7 305 valores, estos datos 
se ingresaron en una base de datos dentro del 
entorno de HEC-RAS. 

Además de los coeficientes de Manning 
correspondientes a cada sección transversal, 
el modelo requiere: a) la separación entre el 
extremo izquierdo de una sección y el extremo 
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Figura 6. Asignación de los valores de Manning según tablas de Chow, Maidment y Mays (1994). 


izquierdo de la sección de abajo; b) la distancia 
que existe entre el cruce de una sección al eje y 
el cruce de la sección de abajo, y c) la separación 
entre el extremo derecho de una sección y el 
extremo derecho de la sección de abajo. Estas 
distancias fueron medidas y se ingresaron en 
una base de datos dentro del entorno de HEC- 
RAS. 


Condiciones de frontera y régimen de flujo 


Una vez capturada la base de datos, se ingresa- 
ron los hidrogramas unitarios del SCS estima- 
dos para cada URH. Dada la sinuosidad de los 
afluentes se eligió simular bajo un régimen de 
flujo mixto, el cual considera las condiciones de 
frontera aguas arriba y abajo del cauce. 


En la confluencia de un afluente tributario 
y el cauce principal no se requirió establecer 
condiciones entre dichas uniones, ya que el 
modelo considera automáticamente condiciones 
de contorno interior, sólo se definieron las con- 
diciones en los extremos aguas arriba y abajo; 
para ello se ingresaron las pendientes a la que 
se asoció un tirante uniforme al inicio y al final 
de cada afluente. 


Índices de evaluación estadística de ajuste del 
modelo 


La calibración busca la correlación entre los 
valores de las variables generadas por el mo- 
delo y las medidas en el sistema real (Ponce, 
1986). Moriasi (2007) recomienda el uso de los 
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siguientes índices estadísticos de prueba para 
determinar la correlación entre ambos valores: 

La raíz del error cuadrático medio (RMSE) 
que se calcula con la siguiente ecuación: 


pl 

N 2 

RMSE = Dala 4) (9) 
N 


La desviación media acumulada (DMA), 
según Molnar (2011), está dada por la siguiente 
ecuación: 


N a,—t,| 


DMA= PIES mt (10) 


Y el coeficiente de correlación (R) dado por: 


q. Fita 
DD aa, a) [Ste Py) 


(11) 


donde a,es el dato estimado por el modelo; t, el 
dato que se asume como medido u observado; 
N, el número de observaciones o estimaciones; 4, 
el promedio de los datos estimados por el mode- 
lo; y £, el promedio de los datos que se asumen 
como medidos (Cai, Liu, Lei, €z Pereira, 2007; 
Cervantes, Arteaga, Vázquez, Ojeda, 6: Queve- 
do, 2013). En los dos primeros índices, el valor O 
corresponde a un ajuste perfecto, mientras que 
un valor más alto representa un menor ajuste. 
Para R, el valor 1.0 indica una alta correlación; 
entre 0.5 y < 1.0, una buena correlación, entre 0.0 
y < 0.5, una baja correlación, y O que no existe 
correlación (ASCE, 1993). 


Resultados 


Al usar el modelo HEC-RAS se simuló el com- 
portamiento hidráulico de los afluentes que 
integran la red hidrológica de la subcuenca 
RH28Be, esto como respuesta a un gasto o escu- 
rrimiento que generaron las lluvias del huracán 
Karl el 17 de septiembre de 2010. 

La simulación hidráulica permitió obtener el 
tirante o nivel del flujo de agua en cada una de 


las secciones transversales, sin embargo, estos 
resultados no pueden ser considerados como 
“aceptables” si no se realiza una prueba de ca- 
libración. Según Pappenberger, Beven, Horritt y 
Blazkova (2005), la importancia de la calibración 
radica en la posibilidad de aplicar sus predic- 
ciones para describir, explicar y comprender un 
sistema real. 

Para llevar a cabo este proceso, se efectuó 
antes un análisis de sensibilidad que consistió 
en realizar diferentes simulaciones variando 
dos parámetros del modelo: a) distanciamiento 
entre secciones transversales y b) coeficiente de 
rugosidad. 

Los resultados obtenidos por las diferentes 
simulaciones (en específico niveles o tirantes) 
se compararon con los niveles observados en la 
estación hidrométrica 28003 (entre las 00:00 y las 
23:00 horas del día 18 de septiembre de 2010); 
para ello se utilizaron tres índices estadísticos 
de prueba que determinaron la correlación entre 
ambos valores. 

La corrida o simulación hidráulica elegida 
para los fines del presente estudio fue aquella 
en la que se obtuvo el mejor ajuste —al aplicar 
los tres índices de evaluación estadística— entre 
valores observados versus simulados (cuadro 3). 

Con la información del cuadro 3 se ajustó un 
modelo lineal entre los niveles medidos y obser- 
vados que se pesentan en la figura 7, donde hay 
una menor relación cuando el nivel supera los 
11 m;, el coeficiente de determinación (R?) para 
el conjunto de datos es de 0.896. 

Las técnicas gráficas de evaluación son herra- 
mientas que facilitan la comparación visual de 
los datos simulados versus observados y permi- 
ten una valoración inmediata de los resultados 
(ASCE, 1993); sin embargo, los métodos gráficos 
revelan poca información sobre la calidad de los 
resultados y el funcionamiento del modelo. Para 
evaluar la capacidad de HEC-RAS se emplearon 
tres índices estadísticos de prueba: RMSE, DMA 
y R. El cuadro 4 muestra el resultado de estos 
índices utilizados en la prueba de calibración. 

Los valores obtenidos indican una buena 
correspondencia entre datos simulados y 
observados, esto confirma que la modelación 
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Cuadro 3. Niveles observados versus simulados (18 de septiembre de 2010). 


Nivel Nivel Nivel Nivel Nivel Nivel 
Hora E e iS Hora observado simulado Hora observado simulado 
observado (m) | simulado (m) 
(m) (m) (m) (m) 
00:00 11.50 13.38 08:00 12.70 11.32 16:00 9.70 10.18 
01:00 13.20 14.40 09:00 10.56 10.85 17:00 9.40 9.92 
02:00 14.80 15.08 10:00 10.10 11.96 18:00 9.16 9.68 
03:00 16.10 15.15 11:00 11.30 11.79 19:00 8.96 9.46 
04.00 15.80 14.46 12:00 11.08 11.54 20:00 8.75 9.32 
05:00 14.90 14.13 13:00 10.78 11.09 21:00 8.60 9.13 
06:00 14.38 13.51 14:00 10.30 10.80 22:00 8.44 9.09 
07:00 13.36 13.02 15:00 10.00 10.46 23:00 8.40 9.09 
16 
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Figura 7. Niveles observados en la estación hidrométrica 28003 Cardel y simulados por HEC-RAS. 


Cuadro 4. Resultados del análisis de comparación de niveles observados y simulados por el modelo HEC-RAS. 
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Índice de evaluación estadística Valor 
Raíz del error cuadrático medio 0.866 
Desviación media acumulada 0.744 
Coeficiente de correlación 0.946 
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hidráulica con HEC-RAS es aplicable para de- 
limitar las zonas con riesgo a inundación de la 
subcuenca RH28Be. 


Delimitación de las zonas con riesgo a 
inundación 


Quaas (2004) menciona que el riesgo depende de 
dos factores fundamentales: el peligro y la vul- 
nerabilidad. Guevara, Quaas y Fernández (2004) 
definen al peligro como la probabilidad de que 
ocurra un fenómeno potencialmente dañino (o 
factor perturbador) durante un cierto periodo 
de tiempo y en un sitio dado. Por otra parte, 
Cardona (2001) define vulnerabilidad como 
la susceptibilidad o propensión de un sistema 
expuesto a un fenómeno o factor perturbador. 
El Cenapred (2004) define al factor perturbador 
como un fenómeno natural o generado por el 
hombre que ocasiona daños a un sistema en un 
grado tal que constituye un desastre. 


Establecimiento del factor perturbador 


Para determinar el peligro se estableció como 
factor perturbador el área indudable (o llanura 
de inundación) obtenida a través de la simu- 
lación con HEC-RAS. En la figura 8 se observa 
que el factor perturbador no representa un 
peligro para todas las comunidades, esto se 
debe principalmente al relieve de la subcuen- 
ca y la distancia que existe entre el centro de 
población y el río La Antigua o afluentes. La 
probabilidad de que ocurra un fenómeno o 
evento potencialmente dañino para una po- 
blación no depende sólo del comportamiento 
hidráulico del cauce. 

De las 96 comunidades que existen en la 
subcuenca RH28Be (donde habitan 34 877 per- 
sonas, según el Censo de Población y Vivienda 
2010), sólo para 23 de ellas representa un peli- 
gro el desbordamiento del río La Antigua (ver 
cuadro 5), fenómeno representado por el nivel 
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Figura 8. Factor perturbador y peligro en un tramo del río La Antigua. 
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O tirante que el agua alcanzó en las secciones 
trasversales durante la simulación hidráulica. 

La frecuencia o recurrencia de inundaciones 
depende de la frecuencia de precipitaciones 
excepcionalmente fuertes; un estudio realiza- 
do por la Conagua (2011) indica que el gasto 
generado por las lluvias del huracán Karl tiene 
un periodo de retorno de 500 años, por lo cual, 
según Ahmed (2008), su frecuencia sitúa a la 
mayoría de las comunidades (de RH28Be) en un 
nivel bajo de amenaza. Esta forma de evaluación 
se usa por lo común en países como Australia, 
Estados Unidos y Suiza; en el presente trabajo se 
utilizó la metodología del Cenapred (2011), que 
aplica como criterios fundamentales el peligro 
y la vulnerabilidad. 


Determinación del índice de vulnerabilidad 


Según Cardona (2001), la vulnerabilidad fí- 
sica corresponde al análisis de las viviendas 
para conocer qué tan susceptibles son a una 
inundación. El análisis se hizo con base en la 
tipología de las viviendas en el ámbito local; 
una vivienda es más vulnerable si los mate- 
riales de construcción son de baja resistencia, 
como adobe, cartón o lámina; y será menos 
vulnerable si está hecha de cemento, ladrillo 
o tabique (materiales de mayor resistencia). 
La información referente al tipo de material 
utilizado para su construcción se extrajo del 
Censo de Población y Vivienda 2010, editado 
por el INEGI (2011). 

El índice de vulnerabilidad se determinó só- 
lo para las 23 comunidades a las que representa 
un peligro el desbordamiento del río La Antigua 
y sus afluentes (como resultado de la simula- 
ción hidráulica). De esta forma se obtuvo que 
11 comunidades, donde habitan 988 personas, 
poseen un índice de vulnerabilidad alto y muy 
alto; cinco comunidades, donde habitan 4 737 
personas, tienen un índice de vulnerabilidad 
medio; y siete comunidades, donde habitan 28 
043 personas, poseen un índice de vulnerabili- 
dad bajo y muy bajo. 


Cartografía de riesgo a inundación 


La base conceptual para elaborar el mapa de 
riesgo a inundación se describe con más detalle 
en Cenapred (2011); sin embargo, se exponen a 
continuación los criterios más sobresalientes de 
esta metodología, que derivó en la cartografía 
de riesgo: a) riesgo alto es toda zona cuyo tirante 
es > 1 m y la velocidad es > 1 m/s; b) riesgo 
medio es la zona cuyo tirante está comprendido 
entre 0.4 y 1 m, y la velocidad fluctúa entre 0.4 
y 1 m/s; c) riesgo bajo es la zona cuyo tirante es 
< 0.4 m y la velocidad es < 0.4 m/s. 

Con base en estos criterios y los resultados 
obtenidos en la modelación hidráulica, se esta- 
bleció el nivel de riesgo de las 23 comunidades 
(cuadro 5) y se generó un mapa de riesgo a 
inundación. La figura 9 representa un bosquejo 
de este mapa correspondiente a la Parte Baja del 
río La Antigua (URH IX). El mapa se editó a una 
escala tal que permite observar la distribución 
de las comunidades, el peligro y la vulnerabili- 
dad de éstas a una inundación. 

La modelación permitió demostrar que las 
lluvias asociadas con el huracán Karl provo- 
caron el desbordamiento de los afluentes que 
integran la red hidrológica de RH28Be; sin 
embargo, no todos los sitios afectados repre- 
sentan un riesgo directo para los habitantes de 
las 23 comunidades, debido a que cinco de los 
nueve tramos simulados impactan sólo en zonas 
agrícolas, pastizales, o relictos de bosque y selva 
(los tramos contenidos en las URH II, III, IV, V y 
VII). Este hecho no debe exceptuar la atención 
de dichas zonas, pues según los reportes de 
Valdez y Arreola (2011), en Veracruz —a causa 
del huracán Karl— se perdió casi 100% de las 
cosechas en cultivos como la papaya y el maíz, 
lo que en términos monetarios representó más 
de 100 millones de pesos (mdp), y en el sector 
pecuario se perdieron más de 1 200 cabezas de 
ganado bovino, 11 000 cerdos y 20 000 aves. 

Por otro lado, las URH L VL VIII y IX ad- 
quieren mayor importancia social. En ellas se 
localiza 100% de las localidades con algún nivel 
de riesgo a inundación. Los reportes de Armen- 
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Cuadro 5. Comunidades de la subcuenca RH28Be bajo algún nivel de riesgo. 


Núm. Comunidad Municipio Población total | Latitud | Longitud Nivel de riesgo 
1 José Cardel La Antigua 24 846 19.37 -96.38 Alto 
2 Juan Yobal Zerna Puente Nacional 4 19.29 -96.61 Alto 
3 Los Cuates Puente Nacional 4 19.28 -96.63 Alto 
4 Puente Nacional Puente Nacional 347 19.32 -96.48 Alto 
5 Salmoral La Antigua 753 19.34 -96.35 Alto 
6 Salmoral parte del Sur Paso de Ovejas 1161 19.34 -96.36 Alto 
7 San José de Enmedio Puente Nacional 6 19.30 -96.58 Alto 
8 San José de La Azud Puente Nacional 9 19.30 -96.58 Alto 
9 Yucatán Paso de Ovejas 543 19.34 -96.36 Alto 
10 Colalco Puente Nacional 9 19.33 96.50 Medio 
11 El Coyolar Puente Nacional 390 19.29 -96.61 Medio 
12 El Cuajilote Puente Nacional 258 19.30 -96.58 Medio 
13 Mata Grande Paso de Ovejas 345 19.35 -96.39 Medio 
14 Paso de la Cruz Puente Nacional 7 19.35 -96.39 Medio 
15 San José Chipila Puente Nacional 402 19.29 -96.59 Medio 
16 Santa María Tatetla Jalcomulco IZ 19.27 -96.71 Medio 
17 Santo Niño de Atocha Puente Nacional 5 19.30 -96.66 Medio 
18 Chichicaxtle Puente Nacional 1910 19.34 -96.47 Bajo 
19 El Cedro Puente Nacional 71 19.30 -96.57 Bajo 
20 El Manguito Paso de Ovejas 173 19.36 -96.38 Bajo 
21 Paso Marino Puente Nacional 388 19.33 -96.46 Bajo 
2 Tomás Vallejo Puente Nacional 3 19.29 96.59 Bajo 
23 Xopilapa Tenampa 357 19.28 -96.78 Bajo 
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Figura 9. Mapas de riesgo a inundación correspondiente al tramo Parte Baja (URH IX) de la RH28Be. 
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ta et al. (2011) mencionan que el huracán dejó 
un saldo aproximado de 500 000 damnificados 
y daños estimado en 70 mil mdp. 


Conclusiones 


Los datos de precipitación registrados durante 
el huracán Karl permitieron estimar los hidro- 
gramas de escorrentía superficial asociados 
con dicho evento meteorológico. Asimismo, 
los puntos de verificación en campo de uso de 
suelo y vegetación facilitaron la asignación del 
valor y de Mamning (o coeficiente de rugosidad) 
mediante el uso de tablas con valores preesta- 
blecidos. El MDE de alta resolución permitió 
representar las condiciones morfométricas 
del río La Antigua y afluentes secundarios 
sin requerir un levantamiento topográfico, 
aunque esto último no significa que el trabajo 
de topografía sea innecesario o deba olvidarse. 
Toda esta información procesada para su uso en 
HEC-RAS permitió delimitar zonas de riesgo a 
inundación. 

Los resultados obtenidos por el modelo se 
compararon con datos observados en la estación 
hidrométrica 28003 mediante el uso de índices 
estadísticos de prueba, los cuales indicaron una 
alta correlación entre ambas series de datos, 
respaldando de esta forma la aplicabilidad del 
modelo HEC-RAS para predecir la magnitud de 
una inundación, catástrofe muy frecuente en la 
subcuenca RH28Be. 

Según los resultados obtenidos en la mo- 
delación hidráulica del río La Antigua y sus 
afluentes, así como el uso de la metodología 
del Cenapred (2011) para elaborar mapas de 
riesgo a inundación (que aplica como criterios 
fundamentales el peligro y la vulnerabilidad), 
se obtuvo que nueve localidades, donde habitan 
27 673 personas, poseen un nivel de riesgo alto; 
ocho un nivel medio con 3 193 personas, y seis 
en un nivel bajo con 2 902 personas. Estos datos 
indican que 96.8% de la población total que 
habita en la subcuenca RH28Be radica en una 
zona con riesgo a inundación. 

De la elaboración y análisis del mapa de 
riesgo a inundación se concluye que las comu- 


nidades más cercanas al río La Antigua y cuyas 
viviendas en su mayoría son de materiales de 
baja resistencia, son las de mayor riesgo. Por el 
contrario, las localidades más alejadas al río La 
Antigua y cuyas viviendas son de materiales 
resistentes, se clasifican en localidades de bajo 
riesgo. 

Cualquiera que sea el caso, desbordamientos 
en zonas agrícolas o centros de población, se re- 
comienda la implementación de medidas estruc- 
turales y no estructurales contra inundaciones, a 
fin de prevenir, controlar y / o mitigar los daños 
ocasionados por estos eventos meteorológicos 
extremos. 
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